Zeeman splitting of zero-bias anomaly in Luttinger liquids 
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Tunnelling density of states (DoS) in Luttinger liquid has a dip at zero energy, commonly known 
as the zero-bias anomaly (ZBA). In the presence of a magnetic field, in addition to the zero-bias 
anomaly, the DoS develops two peaks separated from the origin by the Zeeman energy. We find 
the shape of these peaks at arbitrary strength of the electron-electron interaction. The developed 
theory is applicable to various kinds of quantum wires, including carbon nanotubes. 
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Interaction between electrons in a conductor leads to 
the formation of an anomaly in the tunnelling density of 
states (DoS) at the Fermi level. The one-particle DoS is 
directly related to the differential conductance of a tunnel 
junction, and the anomaly in DoS translates to the zero- 
bias anomaly (ZBA) of the tunneling conductance. This 
anomaly gets stronger if the conductor is disordered, and 
if the dimensionality of the electron system is reduced. 
The perturbative treatment of the DoS anomaly in dis- 
ordered conductors is well-developed . In a disordered 
wire or film, the perturbation theory in the interaction 
strength is divergent at the Fermi level, and therefore 
a non-perturbative treatment is needed to describe the 
DoS at low energies (2[|^,Q]. In one-dimensional conduc- 
tors with one or a few propagating electron modes the 
suppression of the DoS due to the electron-electron repul- 
sion is strong even in the absence of disorder. The den- 
sity of states in this case is adequately described within 
the Luttinger liquid theory ||. The ZBA was observed 
in experiments with higher-dimensional disordered sys- 
tems HQ]. The recently measured [|]|]] strong suppres- 
sion of the tunnelling in a single- wall carbon nanotube at 
low bias gave an evidence that electrons in a nanotube 
indeed form a Luttinger liquid. 

Zero bias anomaly thus provides important informa- 
tion about strongly correlated electron systems. How- 
ever, ZBA is sensitive mostly to the dynamics of electron 
charge but not to that of spin. To probe the spin physics, 
one may study the effect of a magnetic field on the prop- 
erties of the electron system. The perturbative calcula- 
tion shows 0] that the application of a magnetic field 
modifies the anomaly in the DoS. It acquires, in addition 
to the zero-bias dip, two peaks at energies e = dzg^sB, 
where gfJ-BB is the Zeeman energy. The peaks heights 
are equal and proportional to the electron-electron inter- 
action constant in the triplet channel [jjj , which is not ac- 
cessible in a measurement of the conventional ZBA. The 
described Zeeman splitting of the ZBA in disordered nor- 
mal conductors was not observed yet, but a related phe- 
nomenon was studied in thin superconducting films pOl . 



The perturbative theory of Zeeman splitting of ZBA 
in a clean ID system was developed in [0. There, 
the anomalies in the conductance were ascribed to the 
physics of Bragg reflection of electrons off the Friedel os- 
cillation. Magnetic field B splits the standard Friedel 
oscillation in two. The difference between the two cor- 
responding wave vectors is proportional to B. Electron 
scattering off the two components of the Friedel oscil- 
lation results in the conventional DoS anomaly at zero 
energy and two additional peaks in the DoS at ±gp,sB. 
This single-electron picture is valid for a weak electron- 
electron interaction only, and is not applicable to the Lut- 
tinger liquid with a strong repulsion between electrons. 
However, the strongest manifestation of the Luttinger 
liquid behavior is found in carbon nanotubes, where the 
interaction is not weak. Thus, one may question the ex- 
istence of peaks in the DoS at Zeeman energy in such 
systems. 

In this paper, we demonstrate that the tunneling den- 
sity of states in a Luttinger liquid is singular at energies 
e = ±g* [IbB. The effective Lande factor g* here is renor- 
malized by the interaction. The overall magnitude of the 
singular correction to DoS is proportional to the con- 
stant of electron-electron interaction in the triplet chan- 
nel. The energy dependence of the DoS around the sin- 
gularities is given by a power law, 5i/(e) ~ |e ± g*^BB\' y . 
We calculate the exponent 7 in terms of the Luttinger 
liquid parameters. 

Consider the one dimensional Luttinger liquid filling 
the half-line x > and confined by a barrier at a; = 0. 
We decompose the electron creation operator (x) into 
the left- and right- moving parts: = vb\_ s +ip_ s . Here 
± denote left and right movers, and s — ±1 denote two 
spin states. Then, we bosonize electrons p2[: 
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Bosonic variables ifi and ^ describe charge (i = p) and 
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spin (i = a) fluctuations, and a is the short-distance 
cutoff. The fields cpi and i?j are canonically conjugate: 
[(fii(x),-dj(x')] — iSij9(x — x'). Since the electron wave 
function is zero at the barrier [x = 0), the fields <pi sat- 
isfy boundary condition: 



tp p (0) = <p a (0) = . 
The Hamiltonian can be divided into four parts: 



H = H n 



■H' + Hr 
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where the first two terms include the density-density in- 
teractions: 



dx 



(4) 



and H' represents spin-flip backscattering: 

E> = f dx COS\/8</? CT . 
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The last term in Eq. (||) describes Zeeman splitting: 

gH B B f gp B B f 

H B = — J dxp spin (x) = - ^_ J dxd x f a , 

where p sp i n (x) is spin density. The parameters Ui, Ki 
and g± can be expressed in terms of interaction poten- 
tial V(x), but here we treat them as phenomenological 
constants. For free electrons, K p = K a = 1, while for 
repulsive interaction K p < 1 and K a > 1. Also, the bare 
parameters K a and g± are not independent. For g±_ -C 1, 
they are related as K a »1 + g±/2iru a . 

For convenience, we absorb the magnetic field term H B 
into quadratic part by shift ip a — > (p a +gp B BxK a /\f2u a . 
Then, the backscattering term transforms into 



^9 f 
H' = j— I dxcos(V8f a 
(27ra) J J 
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with b = 2gfi B BK a -/u a . 

The Hamiltonian (||) decouples into charge and spin 
sectors. While the charge excitations do not interact 
and their Hamiltonian H p is quadratic, the spin sector 
is described by the sine-Gordon model H a + H' . In zero 
magnetic field the constant g± is renormalized at low en- 
ergies 12 , 



9±(D) 



9±(W) 



1 + 



(0) 



where W is the initial, and D is the running band- 
widths, and g±(W) is the "bare" interaction constant. 
The renormalization group (RG) flow occurs along the 
line K a «1 + g±/2irvp towards the fixed point K* = 1, 
g* ± = 0. Finite magnetic field does not affect the RG 
flow for energies larger than gfi B B. For smaller energies, 



K a becomes essentially independent of e, while g± flows 
toward zero [ p^[ . In this way, the non- linear term H' 
is not relevant at large times, and can be treated as a 
perturbation. 

We are interested in the DoS at e — > g/i B B. This al- 
lows us to consider the Hamiltonian H a +H' acting on the 
states within the energy band of the width D somewhat 
exceeding 2gp B B. The constants K a (D) and g±(D) in 
the Hamiltonian (j^) should be renormalized according to 
Eq. (|). 

Before developing a rigorous calculation, we provide a 
hint, where the singularity in the DoS may come from. 
Tunneling of an electron may be viewed as spreading 
of charge and spin densities, which initially at t — 
were formed near the barrier (x — 0). Because of the 
spin-charge separation, the two density perturbations 
propagate independently. The charge propagates freely, 
since the corresponding Hamiltonian H p in is quadratic. 
The propagation of the spin density is affected by the 
backscattering term Eq. (|^). To demonstrate qualita- 
tively its effect, we expand H' in ip a to the second order 
and then derive the linear equation of motion for the 
field ip a . The first-order expansion term only shifts by a 
small amount the solution of that equation. The second- 
order term generates a contribution oc g± cos(bx)ifa- in 
the equation of motion, and leads to the phenomenon of 
Bragg reflection with wave vector 6/2. As the result, the 
backscattered component of ip a oscillates with frequency 
u! z = u a b/2 = K„g[i B B. These oscillations give rise to 
features in the DoS at energies e — ±cu z . 

We start with retarded Green's function 

oo 

G?(x,x , ,e) = -i J dte ist {{^\{x,t) ^ s {x' 
o 

(here {. . .} is anticommutator) and compute the tunnel- 
ing density of states as []l4| 



dxdx' 



G?(x,x',e) . (7) 



For slowly varying tpi(x) and one may neglect their 

derivatives, and differentiate only the factors e ±lkFX in 
the electron operators ([!]). Equation (Q) can be rewrit- 
ten as 



oo 

/(£) = 2^ Re J dt &(*) e<8 * +Sp(-*) e " fa *] 



(8) 



where 

ft(t) = (Texp ■= exp ■= ) (9) 



V2 



V2 



are time-ordered Green's functions of charge and spin, 
and T denotes time-ordering. 
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To compute these correlation functions at g± = 0, we 
express the fields (pi and in terms of bosonic eigen- 
modes a q , a q of the Hamiltonian Because of the 

boundary condition (|J), only odd modes contribute to 



ipi(x,t) = 2J c ?,i sin qx 
9 



jitii qt 



(10) 



iqx ■ 



Here c 9) i = e~ qa / 2 ^nKi/q, and the short-distance cutoff 
a = u,j/D is related to the reduced bandwidth D. The 
summation in Eq. ( [To|) involves wave vectors g > 
where L is the length of the system. One can compute 
the average in Eq. (^|) using the relations 

valid for any operators A and B linear in a q and a q . At 
zero temperature, the only non-zero average is (a q a q ) = 
1, and one finds 



Ui\t\ + ia 



(12) 



Substituting this expression into Eq. (g) and evaluating 
the integral for e -C D, one arrives to the well-known 
formula l|, 

Co , la 



(0)fV> - 



i/^(e) 



r(i + a) 



(13) 



with the anomalous exponent a = (K p 1 + K a x )/2 — 1, 
and the prefactor 



C = - - - 

To develop perturbation theory in g± , we use the stan- 
dard expression |lq| 



Texp 



exp 



-i»M0,0) 



V2 



s 



(S)c 
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where the averaging is performed over the non-perturbed 
Hamiltonian H a . To the first order in g±, the S-matrix 
is 



DC 

5= 1-i J H'(t')dt' . 



Computing the correction 



2ig A 



OO OO 



i(0 ff (t) - ff (O)) 



cos 



(\/80 ff (a;,t / ) + 6a; ) ~ (cos(-v / 8<?!' (7 (a;,t') + ^) 



we use Eqs. ( |10[ ) and (|11|), and arrive to 



<5£ CT (*) 



:v 00 

d£ J dx 




A' 



(27ra) 2 ./ ./ V a 2 + Ax 2 
x (t) [^+ {x,t,t') e* bx + <5£_ (x, t, t') e- lbx 

where 

&g±(x,t,t) = ± 



(15) 



2xu a t 



u a (t — V) =F x + ia sgn(i — t') 
1 



(16) 



it ff i' =F x + ia sgn(t') 



Unlike Eq. (|l2|), the correction 8Q a (t) contains an oscil- 
lating part. We will retain only this part, since we are 
interested in singularities at non-zero energies. The oscil- 
lation originates from the point of enhanced singularity in 
Eq. (|l6|), x = u a t' = u a t/2. The oscillating contribution 
to 5Q a (t) is 



SGajt) _ ig± 
gV»( t ) ~~2u a \u„\t\ 



9(t)e luJ ^ , (17) 



with u) z = g* fisB, and renormalized Lande factor g* — 
K a g. Using Eqs. (||) and (12) at t > D^ 1 , one finds the 
correction to the DoS 
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which is singular at e = ±w 2 . Since Eq. (18) is even in 
s, we consider further only e sa lo z . Integrating over the 
time domain t ~ |e — lo z \~ 1 , we find for the singular part: 



6v(e) 



e- w. 



gJ-(^) 1 
4it a sin 7T7 

r(l + a) fcosf(c\! — 7) for e> 



(19) 



r(l + 7) \cosf(a + 7) e<w z 
with the exponent 7 = a + 2(X cr — 1), i. e. 



2) + 2{K a - 1) 



(20) 



Equation ( |l9[ ) is valid for arbitrarily strong interaction 
in the charge channel, and confirms the existence of sin- 
gularity in DoS centered at energy 

* u 9* , . 9±(D ~ fJ, B B) 
u z =g VbB ; — = H . (21) 

g iTTUa 

At small g±, which implies K a w 1, the main contribu- 
tion to the exponent in Eq. (19) comes from the charge 
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mode. Therefore, the exponents a and 7 of the power- 
law singularities in the DoS at e = and e — uo z respec- 
tively, are nearly identical, 7 w a. The contribution ( |i"9| ) 
was found for zero temperature, and its energy depen- 
dence is non-analytic at any 7. However, it may be eas- 
ily distinguished from the regular part of v{e) only at 
7 < 1. Also, finite temperature T smears the singularity 
at |e — uj z \ ~ T. 

The DoS anomaly ( pjf ) is directly related to the bias 
dependence of the tunneling conductance G(V) between 
a conventional metal and a one-dimensional conductor. 
The corresponding singular contributions are related as 
5G(V)/G = 8v{eV)/v. Tunneling between the ends of 
two identical one-dimensional conductors (such as an in- 
tramolecular junction between carbon nanotubes ||) also 
has a singularity at Zeeman energy. In this case the peak 
at eV — uj z has a different shape, because it is defined 
by the singularities in the DoS both at e = us z and e = 0. 
The conductance can be calculated as G(V) = dl/dV, 
where the tunneling current /(V) between the two con- 
ductors is proportional to the convolution of the two cor- 
responding DoS: 



eV 



I{V) oc / dev{e)v(e- eV) 



(22) 



Calculating this integral, one finds the singular contribu- 
tion to the differential conductance 



SG(V) 
G(V) 



1 



Au a sin I (a 
r(l + 2a) 



•7) 



eV 



a+1 



(23) 
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The exponent a + 7 here coincides with 27 up to a small 
term of the order of K a — 1. This "exponent doubling" 
at eV — ui z is similar to that occuring at zero bias (9) . 

It is interesting to analyse Eq. (23) in the limit of weak 
interactions, in which m2\i 



g± u(2k F ) 



2irvi 



2ttvp 



TTVf 



To the first order in the interaction potential U (q) , 
Eq. (H) yields SG/G = {U{2k F )/^Ku a ) ln(\eV-w z \/u g ). 
This result is in agreement with the first-order expansion 
of the tunneling conductance obtained in [TlJ| . However, 
beyond this order, there is a difference between Eq. j2^ ) 
and Eq. (47) in It stems apparently from the inap- 
plicability of the RG approach developed in |ll|] for the 
treatment of tunneling at energies close to u> z . 

We derived Eq. ( |l9| ) and ( p0| ) for the single-mode Lut- 
tinger liquid. In the case of carbon nanotubes, one has to 
take into account the degeneracy between the two conic 
points in the Brillouin zone Treating the interac- 

tion in the charge channel non-perturbatively, as we did 
before, we find the exponent of the singularity at e = uj z . 



7nt ~ «nt 



(24) 



Thus, the exponent 7 nt again nearly coincides with the 
ZB A exponent a n t . The reported values of a n t in the ex- 
periments with carbon nanotubes were a n t = 0.3 -7- 0.6, 
and therefore the peak at e = ui z should be sharp and 
easy to observe. 

To conclude, the application of a magnetic field to 
a Luttinger liquid creates additional singularities (two 
peaks) in the tunneling density of states. We have 
found the power law characterizing these singularities 
and demonstrated that they are robust, i. e., they persist 
at any interaction strength in the charge channel. The 
magnitude of the peaks is determined by the short-range 
interaction, which plays a minor role in the charge physics 
of a Luttinger liquid and therefore is hardly accessible in 
the measurements of the conventional zero-bias anomaly 
of the tunneling conductance. 
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